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Abstract 

The scaling behavior of pure gauge SU(3) in the region f3 = 5.85—7.60 is examined by 
a Monte Carlo Renormalization Group analysis. The coupling shifts induced by factor 
2 blocking are measured both on 32*^ and 16"^ lattices with high statistics. A systematic 
deviation from naive 2-loop scaling is clearly seen. The mean field and effective coupling 
constant schemes explain part, but not all of the deviation. It can be accounted for by 
a suitable change of coupling constant, including a correction term 0{g'^) in the 2-loop 
lattice /9-function. Based on this improvement, \/a/ is estimated to be 2.2(±0.1) 
from the analysis of the string tension a. 
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Since a confirmation of tlie approacli to tlie continuum limit lias basic importance in 
lattice Quantum Chromo Dynamics (lattice QCD), systematic scaling analyses at larger 
(5 are inevitably required. Recent analyses have shown that scaling violations persist 
in physical quantities such as the string tension and hadron masses up to /? = 6.8.|jl| 
On this problem, it has been argued that a suitably chosen coupling constant reveals 
perturbative scaling. [|, |, |, |, ^ 

Monte Carlo Renormalization Group (MCRG) tells us the coupling shift A/? induced 
by scale transformation of the lattice spacing, a —>■ sa-J 

where fil/g^) = aA and Pf{g) is the lattice P function. (For bare lattice coupling 
constant, l/g"^ = /3/6.) This gives us another way to examine the scaling behavior of 
lattice QCD. In SU(3) lattice gauge systems, MCRG has been performed by several 
groups on 16^ lattices in the large /? region up to 7.2. ^, [1^, |n|, However, these 



results were not conclusive and even controversial. Although Gupta et al.|]rT| have 
claimed consistency with asymptotic scaling from their a/3 blocking result. Bowler 
et al.p, have found sizeable deviations from 2-loop scaling. Recently Hoek has 
reanalyzed the same data and claimed a very slow approach to scaling. [|12| These varied 
conclusions are mainly due to the difficulty of obtaining a precise value for the coupling 
shift in high (3 region, because the deconfining transition prevents an accurate matching 
between Wilson loops on a small lattice. Studies on larger lattices with better statistics 
are required to clarify the scaling behavior in the high /3 region. 

In this work, we report results of an MCRG study and scaling analysis in the high 
/3 region both on 32^* and 16^ lattices with high statistics. 

Using a 32^ lattice gives us the following advantages. 

i) one more blocking level than previous works. It is not apriori clear that blocking 
from 16*^ to 2^ is sufficient for loop matching. A deeper blocking is preferred to confirm 
good matching and for closer matching conditions. 

ii) we remain in the confinement phase up to /? ~ 6.9 whereas a 16^ lattice is above 
the deconfining transition point for (3 > 6.35. 



The numerical simulation has been performed on 512-cell parallel processor Fujitsu 
AP1OOO.[0 Three lattices of size 32^, 16^ and 8^ are generated by the over-relaxed 
pseudo heatbath algorithm IT^ (the mixing ratio of over-relaxation to pseudo heat bath 
is 9:1 in average). Blocking is performed every 10 updates. 2K - 3K configurations are 
blocked at each (3 (more near the deconfining transition point). For error estimation, 
the jack-knife method is applied. We monitor autocorrelations of blocked Wilson loops 
in those measurements to know the statistical validity of sampling. (See Fig.3) Near the 
deconfining transition point, twice more configurations are blocked. Effect of long au- 
tocorrelation there is taken into account in the error estimation. Details of prescription 
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will be presented elsewhere. 

In this work, the Swendsen blocking transformation is used to double the lattice 
spacing(s = 2)|]T5|. Blocking is repeatedly performed down to 2"^ for two lattices, one 
of size L { aX (3) and the other of size L/2 (at f5 — AjS). To match long range physical 
contents on both lattices, a set of Wilson loops on one blocked lattice, is compared 
with the corresponding one on the other blocked lattice. For early matching between 
blocking trajectories, the blocking transformation is controled by a parameter q which 
governs the size of Gaussian fluctuations around the maximal SU(3) projection of the 
block link variable ||16[. The coupling shift Af3 is determined at the value of q where 



the mismatch of the two sets of Wilson loops is minimum. Planar 1x1 and 1x2 
Wilson loops, and non-planar 6-link ("twist" and "chair") and 8-link ("sofa") loops are 
measured on blocked lattices at g = 0.0,0.02,0.04 for this purpose (see Fig.l). A brief 



description is found in our previous report] 13 



AjS a.t (3 = 5.85 — 7.60 measured by matching 16^ and 8^ lattices after three blocking 
steps is shown by open circles in Fig. 2a. Improved statistics give a clear systematic 
behavior for the new data in comparison with those of previous works (Fig. 2b) although 
they are consistent. The data are significantly below the 2- loop scaling result (solid 
curve). Therefore naive 2-loop scaling does not hold in this region. The deviation 
rapidly decreases as (3 increases. But, even at /5 = 7.60, 10% deviation remains. 

A notable feature of the present data is the following. For (3 < 6.3 , the quality 
of the data is sufficient and shows the approach to the 2-loop scaling result. On the 
other hand, the data above /3 = 6.3 suffer relatively large errors. The matching of 
blocked Wilson loops becomes difficult in this region. This difficulty is caused by the 
deconfining transition on the 16^ lattice. As shown in Fig. 3, the autocorrelation time 
of blocked Wilson loops sharply peaks at /3 = 6.35 ±0.05 ( On 8"^ lattice, /3 = 5.90). An 
increase of fluctuations of the blocked Wilson loops is also evident at this point. Thus, 
the 16^ lattice turns into deconfining region at this point. It is noted that the value 
(3 = 6.35 is smaller than the previously expected value 6.45 ± 0.05 from an analysis of 
the finite temperature phase transition |18[. Above this f3, Wilson loops are dominated 
by perturbative contributions and the matching suffers large errors. 



Matching between 32^ and 16^ has been tried at (3 = 6.35 — 7.00. The deconfining 
transition point on a 32^ lattice is pushed up to around P = 6.9 since the coupling 
shift is ~ 0.54 in this region and the deconfining transition point is (3 = 6.35 on a 16^ 
lattice. At /3=6.35,6.55,6.65 and 6.80 where the lattice is in the confinement phase, 
the matching can be performed successfully and the resultant coupling shifts after 3 
and 4 blockings agree within error as shown in Fig. 4. Thus,the data become stable 
for blockings greater than three and this fact gives reliability to the 16^ lattice (three 
blockings) measurements. 

The measured coupling shifts (black circles in Fig. 2a) are naturally connected with 
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those of (3 < 6.35 on the 16^ lattice. Thus, a systematic deviation from naive 2-loop 
scahng in this region is further confirmed by these measurements. 

At (3 = 7.00 where the 32^ lattice is above the deconfining transition point, definite 
data could not be extracted from our 3K blocked configurations (30K sweeps) due to 
very long range fluctuations. 



Although a trend to approach the 2-loop scaling value is seen, the coupling shift 
shows significant deviation. This deviation can be partly absorbed by a mean field 
scheme ^ or an effective coupling constant scheme. 0, ^ |^ In the mean field scheme, 
the MS coupling constant is given by l/gjjgij^jaY =< Upiaq > / g'^ + 0.025 where 
< Upiaq > is the average plaquette. Similarly the effective coupling constant is defined 
as (yfg = 3(1— < Upiaq >)• Both coupling constants are obtained here based on our 
measurement of < Upiaq > at /? = 5.70 — 7.60. The coupling shift is calculated assuming 
the 2-loop scaling form for ^xTsP^ 



/2(x)=vr(— A^)^exp(-— ), x = ^ (2) 

where bo = 33/(487?^) and h = (102/121)6^. As shown by curve M in Fig. 5, this scheme 
partly explains the present MCRG data. Similary, the effective coupling constant ge 
improves the agreement somewhat, as shown by curve E in Fig. 5. 

However, a sizeable deviation still persists between these schemes and the present 
data. Therefore we need a different prescription to approach the continuum limit from 
the presently accessible region of (3. In ref.2,3, to get the continuum limit, a linear 
or bilinear extrapolation in terms of the lattice constant a was assumed for physical 
quantities expressed in unit of A^g^, while the authors of ref.6 used a linear extrapolation 
in 1/lna. Here instead, the continuum limit is extracted by taking into account the 
lattice 13 function and expressing our MCRG data in terms of an effective coupling 
constant. We fit the data by the lattice /3-function including a next-order correction, 
and define our effective coupling l/gf^ by a free shift of l/fi'lj;^ as follows: 

dxu , bi b' , . 

&o + — + ^, (3a) 



with 

Xu = — = ^— - Xq . (4) 

Assuming the correction term b' is small, we actually use the following equation instead 
of eq.(3a), 

dxu 1 b' 
7—\-, 1 — r — — n] = dlna. (3b) 

The solution of eq.(3b) is 
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f{xu) = /2(a;«)exp(— 2— ) . (5) 

Using eqs.(4) and (5), we have attempted two fits for the data above /3 = 6.00 : (A) 
Restrict Xu = l/gjjg, (i.e.Xo = 0) with b' as a free parameter, (B) Allow a shift of 
the effective coupling constant, i.e. both xq and b' are free parameters. The results 
are shown in Fig.5. In case (A), only a poor fit is obtained. The coefficient b' is also 
relatively large as b' /bo = -0.0850(0.0013) while 61/60 = 0.0587. Thus, it is difficuh 
to explain the present data by an additional 0{gjj^) term only. On the other hand in 
case (B), we can fit the data quite well by values of parameters as (curve B) 

b' 

xo = 0.442(0.004), — = -0.0119(0.0008) . (6) 

bo 

In this case, the coefficient of the 0{gD term is small but the shift parameter xo is 
non- negligible. It is noted that we have tried fitting in the effective coupling scheme 
also and get similar results. Since the sign of b' is negative, we say that our lattice is 
coarser than predicted by 2-loop scaling. 

Based on this scheme, we can discuss the continuum limit of physical quantities in 
units of A-g^ . In Fig. 6, scaling of the string tension of ref.s 6,20,21 is examined. As 

expected, good scaling is obtained for ^ / iC^^ for /3=5.7-6.8 as shown by black circles 
in the figure and the value is 

V^/A^° = 2.2(0.1). (7) 

This value is also interpreted as the continuum limit value in the present scheme. Our 
result (7) is larger than that of the static quark potential in ref.6 (1.80 ±0.06) and ref.22 
(1.72 ± 0.13). We note also that A^° = 200Me1/ for ^ = UOMeV. It is slightly 
smaller than that extracted from Ip — Is splitting of charmonium in ref.2(234Mey). 



In this work, the scaling behavior of SU(3) lattice theory in the interval (3 = 5.85 — 
7.60 is studied by MCRG analysis. A significant deviation from naive 2-loop scaling 
is seen. This is clearly observed in the data by measurements kept in the confinement 
phase up to /5 = 6.80. Thus naive 2-loop scaling does not hold in this region. It is 
shown that a large part of the deviation is accounted for by the mean field and/or 
the effective coupling constant schemes. We show further that the deviation can be 
absorbed by a next-order correction to the 2-loop /3-function, together with a shifted 
mean field coupling constant. The coefficient of the next order term, 6', is consistently 
small. Based on the scaling behavior of this lattice /5-function, the string tension remains 
constant in the region (3 = 5.70 — 6.80 where the lattice spacing changes by a factor 
six. The estimated A-^ in the continuum limit is 0.46(0.03)a/o". Though there is 
still a 20% discrepancy between different methods^, y, this seems an encouraging 
confirmation of the approach to the continuum limit. It is also noted that the suggested 
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coupling constant Qu diverges at some value of the bare coupling constant( [3 ~ 4.5). An 
interesting possibility is that this divergence is related with the transition point from 
strong to weak coupling region studied by Bhanot and Creutz in a space of couplings 
of Wilson action and that of adjoint representation. [p3 



Finally the significance of MCRG measurements in confinement phase is stressed. 
Above the deconfining transition point, the matching procedure inevitably suffers large 
errors. Although most of the data measured above the transition are consistent, those 
in confinement phase have higher quality and reliability. Study on larger lattice is 
apparently preferred. In order to cover the region up to (3 =7.5, we need a 64'* lattice. 
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Figure Captions 

Fig.l Wison loops used in lattice matching. 

Fig.2 Coupling shift A/3 measured in this work(a) and that in previous works(b). Solid 
curve shows 2- loop asymptotic scaling with bare coupling constant. 

Fig.3 Autocorrelation time r on 16^ and 32^ lattices. 

Fig. 4 Coupling shift Aj3 on 32^ lattice in different blocking level. 

Fig. 5 Coupling shift reproduced by different coupling constant schemes. Notations are 
described in the text. 

Fig.6 Scahng of string tension in different coupling constant schemes. 
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